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ABSTRACT 


The  classic  single  echelon  inventory  model  is  restruc¬ 
tured  as  a  two  echelon  problem  in  which  demand  and  resupply 
are  deterministic.  Using  cost  minimization  as  the  objective, 
three  models  are  developed  which  address  the  problems  of 
(1)  no  stockouts  allowed  (EOQ) ,  (2)  backorders  allowed,  and 
(3)  finite  production  with  no  stockouts  allowed.  General 
solutions  for  the  optimal  policy  are  obtained  in  the  EOQ  and 
finite  production  models.  In  the  backorder  model,  the  ana¬ 
lytical  argument  is  limited  to  the  case  in  which  only  time 
dependent  backorder  costs  occur.  Algorithms  are  developed 
for  solving  problems  for  all  three  models,  and  selected  para¬ 
meter  values  are  used  to  test  the  behavior  of  the  models. 
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I.  INTRODUCTION 


Traditional  inventory  analysis  addresses  the  single 
echelon  inventory  problem.  In  this  problem,  customer  demands 
are  received  and  issues  are  made  from  a  single  outlet.  The 
outlet  in  turn  is  replenished  from  a  single  source. 

Unfortunately,  most  real  world  inventory  systems  are  not 
this  simple.  A  major  manufacturer,  for  example,  usually  has 
a  production  and  distribution  system  that  includes  (1)  fac¬ 
tories,  (2)  factory  warehouses,  (3)  regional  warehouses, 

(4)  local  warehouses,  and  (5)  retail  outlets  [Hadley  and 
Whitin] .  Such  a  system  is  an  example  of  a  multiechelon  inven¬ 
tory  system.  Each  level  of  the  system  is  a  separate,  distinct 
echelon.  Figure  1  depicts  graphically  this  organization.  The 
United  States  Navy  with  its  organization  of  inventory  control 
points,  supply  centers,  supply  depots,  shop  stores,  tenders  and 
underway  replenishment  ships  provides  an  example  of  a  large 
multiechelon  inventory  system. 

In  defending  their  single  echelon  models  Hadley  and  Whitin 
argue  that,  even  though  most  real  world  systems  are  multi¬ 
echelon,  it  is  often  true  that  the  system  need  not  be  treated 
as  multiechelon.  They  contend  that  a  different  organization 
frequently  operates  each  level  (echelon)  of  the  system.  In 
the  example  of  Figure  1,  a  factory  and  its  associated  ware¬ 
houses  might  be  operated  by  an  equipment  manufacturer  whereas 
the  regional  warehouses  might  be  separate  privately  owned 
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Figure  1.  Typical  Multiechelon  Inventory  System 

distributorships  not  under  the  administrative  control  of  the 
manufacturer.  In  like  manner#  the  local  outlets  might  be 
individually  owned  retail  stores.  Thus  we  would  have  a  sys¬ 
tem  comprised  of  several  single  echelons  which  are  linked  by 
a  physical  dependence  but  separated  by  administrative 
controls. 

While  it  is  clear  that  there  are  many  examples  of  multi¬ 
echelon  systems  that  can  be  treated  as  a  series  of  independent 
single  echelon  problems,  it  should  be  equally  clear  that  there 
is  a  large  class  of  systems  which  must  be  treated  as  strongly 
dependent  multiechelon. 

It  is  the  objective  of  this  thesis  to  develop  and  investi¬ 
gate  the  behavior  of  multiechelon  models  constructed  from  the 
classic  single  echelon  inventory  models.  Specifically,  the 
following  deterministic  single-echelon  models  of  Hadley  and 
Whitin  will  be  restructured  as  two-echelon  models: 
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II.  THE  TWO  ECHELON  MODEL  WITH  NO  STOCKOUTS 

In  any  study  of  inventory  theory,  one  of  the  first  models 
investigated  is  the  deterministic  demand  model  with  no  stock¬ 
outs  allowed;  often  referred  to  as  the  Economic  Order  Quantity 
(EOQ)  model.  While  it  is  true  that  complete  deterministic 
demand  is  almost  never  known,  it  is  felt  that  the  mathematics 
of  this  model  provide  a  good  starting  point  for  any  inventory 
analysis.  Further,  it  is  felt  that  the  deterministic  model 
will  provide  an  insight  into  the  operation  of  a  stochastic 
demand  model.  It  is  considered  appropriate,  therefore,  to 
begin  this  multiechelon  analysis  by  considering  a  two  echelon 
extension  of  the  classic  EOQ  model. 

A.  MODEL  FORMULATION 

The  two  echelon  model  is  based  on  the  following 
assumptions : 

1.  The  upper  echelon  always  replenishes  its  supply  from 
the  same  outside  source.  The  lower  echelon  always  replenishes 
its  supply  from  the  upper  echelon. 

2.  The  upper  echelon  is  always  able  to  meet  the  demand 
of  the  lower  echelon, 

3.  External  customer  demand  always  occurs  at  the  lower 
echelon. 

4.  The  external  demand  rate  is  deterministic,  continuous, 
and  constant  with  a  value  of  X  units  per  year. 
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5.  Neither  backorders  nor  lost  sales  are  allowed. 

6.  Procurement  lead  time  (PLT)  is  negligible. 

7.  The  order  quantity  at  the  upper  echelon,  Q2,  is  an 

9  f 

integer  multiple  of  the  lower  echelon  order  quantity,  Ql. 

Let  n  =  Qj/Q-^  '  ^en  this  assumption  requires  n  >  l  and 
integer. 

The  optimal  policy  throughout  this  thesis  will  be  that 
which  minimizes  average  annual  variable  system  cost,  subject 
to  the  constraint (s)  of  the  model.  The  form  of  the  average 
annual  variable  system  cost  will  be  developed  by  first  deter- 

9 

mining  the  system  cost  per  cycle,  and  then  dividing  this  cost 
by  the  cycle  length.  The  total  system  variable  cost  per  cycle 
is  the  sura  of  the  individual  echelon  variable  cycle  costs, 
K^(c)  and  K2(c).  The  cycle  length  is  defined  as  the  time  be¬ 
tween  receipt  of  two  successive  orders  at  the  upper  echelon 
(echelon  two) . 

The  cost  per  cycle  at  echelon  two  is  comprised  of  an 
ordering  cost,  an  inventory  holding  cost  (IHC)  and  a  purchase 
cost.  The  ordering  cost  is  assumed  to  be  a  constant  cost  per 
order  which  includes  the  administrative  costs  associated  with 
inventory  review  and  order  (contract)  preparation.  This  cost 
is  independent  of  the  quantity  on  hand  or  on  order.  The  IHC 
includes  a  warehousing  cost,  an  obsolescence  cost  and  a  fore¬ 
gone  opportunity  cost.  It  is  assumed  to  be  a  function  of  the 
inventory  on  hand  and  to  be  expressable  by 

IHC  =  l,c  /c£,(t)dt, 

0  c 


On  hand  inventory 


~,^T?yy$yv;v  *&■  ;  ^  ?  y ;  ps^c 


T 


k2(c)  =  A2  +  CQ2  +  I2c  u  2  (t)dt 


(2-1) 


The  height  of  the  on-hand  inventory  at  echelon  two  is 

Q2  -  Qi  =  tn-1]  Qx 

at  the  beginning  of  the  cycle  since  Q2  is  ordered,  arrives 
immediately,  and  is  partially  used  to  meet  an  immediate  demand 


of  Q1  units.  The  next  demand  at  echelon  two  is  for  an  amount 


which  occurs  Q^/X  time  units  later.  Similarly,  another 
demand  of  occurs  after  Q^/A  additional  time  units  have 
passed.  The  on-hand  level  at  echelon  two  continues  to  de¬ 
crease  in  steps  of  Q1  until  it  reaches  zero.  It  remains  at 
zero  for  the  last  Q^/x  time  units  of  T.  The  area  under  the 
on-hand  inventory  curve  for  one  cycle  of  echelon  two  is 
therefore 
.t 


J^l2(t)dt  =  r-J/Oi  »  Qi  • 


From  (2-1) ,  the  total  variable  cost  per  cycle  at  echelon 
two  is 


K2(c)  =  a2  +  cq2  + 


I2C  n [n-1] 


2X 


(2-  2) 


From  the  assumption  that  Q2  =  nQ^  ,  where  n  is  integer, 
it  follows  that  the  cycle  length  at  echelon  two  is  equal  to  n 
reorders,  at  echelon  one.  Therefore,  from  Figure  2, 


nQ, 


T 


JT 


‘1 
T~  * 


(2-3) 


The  average  annual  cost  at  echelon  two  is  obtained  by 
dividing  (2-2)  by  (2-3)  which  yields 


12 


■A 
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A-A  I-C[n-l]Qf 

2  +  cx  +  2 - L 


or,  as  a  function  of  and  n, 

A2A  I2CQ1  [n-1] 


+  CA  + 


(2-4) 


Development  of  the  cost  per  cycle  at  echelon  one  follows 
an  argument  that  is  analogous  to  the  echelon  two  development, 
except  that  there  is  no  purchase  cost  incurred  at  echelon  one. 
Also,  from  the  assumption  that  Q2  is  an  integer  multiple  of 
Qx  ,  it  follows  that  in  each  cycle  there  are  n  reorders  at 
echelon  one,  i.e.,  K^c)  is  linear  in  n.  Thus  the  cost  per 


cycle  can  be  written  as 


K^fc)  =  nA^  + 


nIi®  l 


(2-5) 


When  (2-5)  is  divided  by  (2-3)  the  average  annual  variable 
cost  is 


A1X  11CQ1 
qt  +  — 


(2-6) 


The  total  average  annual  cost  of  the  system,  K,  is  the  sum 


of  (2-4)  and  (2-6). 

A,A  I.CQ-,  A0A  I0C  [n-1]  Q, 

K  -  oj-  +  -V1  +  +  cx  +  -2-s - 1 


(2-7) 


The  model  has  assumed  that  -he  unit  price  is  independent 
of  the  quantity  ordered.  Therefore  the  C  term  in  (2-7)  is  a 
constant,  and  can  be  dropped  from  (2-7)  resulting  in  (2-8) 
being  the  average  annual  variable  cost. 
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AXX  IjCQj.  A2X  i  2C  [n-1]  Q^ 


(2-8) 


B.  OPTIMAL  POLICY 

Determining  the  optimal  inventory  policy  involves  finding 
the  values  of  and  n  (call  them  Q*  and  n*)  which  minimize 
(2-8)  for  a  given  set  of  model  parameters.  To  find  this 
policy,  rewrite  (2-8)  as 


A2  X  [I1  +• IoIn-iJlCQ, 

K(Qx.n)  -  t*!  +  ^1  ^  +  -L- j2 - -1 


(2-9) 


Then  since  K(Q^,n)  is  convex  in  Q1  for  fixed  n,  Q£(n)  must 
satisfy  the  equation 

3K  „  .  A2iX  .  tJl  +  I2ln“1JJC 

^  =  “  tAl  +  JT3^2  +  - 2 -  =  0  ' 


Qj;(n)  = 


A-  “11/2 

2X  tAl  +  IT3 


“1  T  2 


(2-10) 


After  substituting  (2-10)  for  in  (2-9)  and  collecting 

terms,  (2-9)  takes  the  following  form: 

A?  1/2 

K (Q*  (n)  ,n)  =  {2X  [Ax  +  ~]  (i!  +  I2(n-l]]C) 


[2XA,I. 
2XA1I1C  + 


+  2XA1I2C[n-l]  +  2XA2I2C[ 


Note  that  the  first  term  in  (2-11)  is  equivalent  to  K' 


(2-11) 


(Eq.  (2-11)  [Hadley  and  Whitin] ) .  In  fact,  it  .is  readily  seen 
that  each  term  in  (2-11)  contains  a  form  of  this  classic 
formula.  Let 

Kij  -  2XAiIjc  • 


Then  (2-11)  can  be  rewritten 


K*  1/2 

K(n)  =  K(Q*  ,  n)  =  {K2  +  —  +  K22[n-1]  +  K22  I2—] } 

(2-12) 

It  is  immediately  seen  that  (2-12)  is  positive  for  all 
feasible  values  of  n  (n  >.  1) .  Further, 

lim  K (n)  =  [K2  +  K2  J1/2  , 

n-*l  w  21 


If  the  assumption  that  n  is  an  integer  is  temporarily  re 
laxed  the  slope  of  K(n)  can  be  readily  investigated  by 


considering 


3K 

3n 


(2-13) 


If  (2-13)  is  positive  when  n  =  1,  then  it  follows  that  it 
is  positive  over  all  n,  and  the  optimal  value  of  n  is  n*  =  1. 

If  (2-13)  is  negative  when  n  =  1,  then  there  exist  values  of 
n  such  that  K(n)  is  less  than  K(l).  However,  since  it  has 
been  shown  that  (2-13)  is  positive  for  very  large  n,  it  follows 
that  n*  is  finite,  and  satisfies 


■'■  ^ ',‘f '5 '■  V^  '  "-*' 


3K 

Tii 


=  0  , 


or 


n’ 


rK2  _  2  ‘ 
*21  *22 


K 


12 


1/2 


(2-14) 


A  further  consideration  of  the  case  where  (2-13)  is  nega¬ 
tive  at  n  =  1  reveals  that 


32K 


3n 


over  the  range  of  n  values  from  1  to  beyond  n*.  Let  the  value 
of  n  where 


fl  =  0 

3n 


be  denoted  by  n.  Then  in  the  range  1  <  n  <  n  the  function 
K(n)  is  convex.  For  values  of  n  >  n,  (2-12)  is  not  convex. 
However,  (2-13)  will  always  be  positive.  Therefore  n*,  given 
by  (2-14) ,  is  the  optimal  non-integer  value  of  n.  The  reader 
will  note  that  the  above  arguments  indicate  that  (2-12)  is 
pseudo- convex . 

Since  K(n)  is  pseudo-convex  over  all  n,  it  follows  that 
K(n)  is  also  pseudo-convex  on  the  integer  values  of  n.  If  n* 
is  not  an  integer,  then  the  optimal  integer  value  of  n  will  be 
either  n  (the  smallest  value  greater  than  n*)  or  n  (the 
largest  value  less  than  n*) .  The  integer  optimal  will  be  the 
value  of  n  corresponding  to  K*  =  min  [K(n)  ,K(n)]  . 

In  summary,  the  optimal  inventory  policy  is  found  as 
follows : 
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1.  Compute  K21,  K^,  K^  . 

2.  Test  K^  -  Kj2  2.  Ki2  *  If  this  inequality,  holds  then 
n  optimal  =1.  Go  to  step  8. 

3.  If  -  Kj2  <  k^2  then  compute  n*  using  (2-14)  . 

4.  If  .n*  is  an  integer  then  n  optimal  =  n*.  Go  to 
step  8. 

5.  If  n*  is  not  an  integer,  compute  n  and  n. 

6.  Compute  K(n)  and  K(n)  from  (2-12). 

7.  K*  =  min  {K(n  ,  K(n)}  and  n  optimal  is  the  value  of  n 
corresponding  to  K*  . 

8.  Compute  Q£  from  (2-10) . 

C.  NUMERICAL  EXAMPLE 

As  an  illustration  of  the  algorithm,  consider  the  follow¬ 
ing  problem  where 


x  =  0.75 

A2  =  $200.00 

2  =  0.50 

\  -  100  units/year 

^  =  $25.00 

C  =  $100/unit 

.  .  2 
Begin  by  computing  =  2XA^IX  • 


k|x  =  3.0  x  106 
=  2.0  x  106 
K^2  =  2.5  x  105 

Next  compare  K21  -  K22  =  1.0  x  106  >  K^2  =  2.5  x  105  , 
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so  compute  n*  using  (2-14) .  The  value  is  2.0  which  is 
integer.  Therefore  n  optimal  =  2.  Using  (2-10)  the  optimal 
order  quantities  are 

QJ  =  166.67  units 

and 

=  166.67  x  2  =  333.34  units. 

D.  PARAMETER  ANALYSIS 

The  primary  objective  of  this  analysis  was  to  observe  the 
effect  of  parameter  variations  on  optimal  n.  The  values  of 
four  of  the  five  parameters  were  fixed  and  the  value  of  the 
remaining  parameter  was  allowed  to  vary.  The  process  was  then 
repeated  for  each  of  the  other  parameters  so  that  the  effect 
of  each  could  be  observed. 

Representative  graphs  of  optimal  n  versus  a  given  para¬ 
meter  are  illustrated  in  Figures  7  through  10.  In  Figure  7 
it  is  seen  that  n  optimal  is  inversely  related  to  A^.  How¬ 
ever,  as  noted  in  Figure  8,  there  is  a  direct  relationship 
between  n  optimal  and  A2*  Further,  a  comparison  of  Figures 
7  and  8  reveals  that,  for  the  parameters  selected,  n  optimal 
is  much  more  sensitive  to  than  to  A^. 

Figures  9  and  10  indicate  the  relationship  of  n  optimal 
to  I1  and  I2,  respectively.  The  results  indicate  that  n 
optimal  is  directly  related  to  1^,  inversely  related  to  I2/ 
and  that  n  optimal  is  more  sensitive  to  a  change  in  holding 
cost  at  echelon  two  than  at  echelon  one. 
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An  analysis  of  the  effect  of  demand  rate  on  n  optimal 
revealed  that  n  optimal  was  completely  insensitive  to  X  over 
the  range  5  <L  X  <,  10,000. 
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III.  FINITE  PRODUCTION  AT  THE  SECOND  ECHELON 


An  additional  two  echelon  variation  to  the  deterministic 
models  treated  [Hadley  and  Whitin]  addresses  the  problem  of 
finite  production.  This  model  would  have  application  in  any 
situation  where  the  upper  ee'  elon  manufactures  as  well  as 
warehouses  a  given  product. 


A.  MODEL  FORMULATION 

The  model,  shown  graphically  in  Figure  3,  assumes  that 
echelon  two  produces  the  needed  material  at  a  constant  con¬ 
tinuous  rate  of  ip  units  per  unit  time  (ij>  >.  X )  .  The  cost  A2  is 
now  considered  to  be  a  set  up  cost  which  is  incurred  each  time 
a  new  production  run  is  initiated.  This  set  up  cost  is  inde¬ 
pendent  of  the  quantity  produced.  The  rest  of  the  assumptions 
are  identical  to  those  of  the  EOQ  model. 

The  model  will,  as  in  Section  II,  seek  optimal  values  of 
and  Q2  which  minimize  the  average  annual  system  cost.  This 
cost  is  the  sum 

K  =  Kx  +  K2  , 

where  is  the  a^exctge  annual  cost  at  the  ith  echelon.  In 
Section  II  it  ;as  shown  that. 


Ki  ■  -j-- +  ihCQi  ■ 

X 

To  develop  K2  it  is  necessary  to  define  the  cycle  length 


T  =  t.  +  T  +  T,  , 
1  p  d  ' 


(3-1) 
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where  tj,  is  the  time  from  the  beginning  of  the  cycle  until  the 
start  of  production,  T  is  the  production  length,  and  T,  is 
the  time  required  to  deplete  the  inventory  at  echelon  tw> 
after  production  stops  (see  Figure  3) .  Note  that 


T 


Qi 

jr 


w2 

3 r ' 


wi 

v  # 


and 


T. 


Define  i  as  the  number  of  times  that  the  quantity  is 
demanded  at  echelon  two  during  the  production  period,  T  .  It 

tr 

can  be  seen  that 


t,  +  T  ^ 

i  -  t-i-r-E] 


(3-2) 


where  t  = 


Q1 

r* 


As  in  Section  II,  the  variable  cost  per  cycle  at  echelon 
two,  KjfC),  will  have  the  form 


K. 


t,  +  T  +  T, 
1  p  d 


(C)  *  A2  +  I2C  t)  dt 


(3-3) 


From  Figure  4  it  follows  that 


tx  +  Tp  +  Td 


Jj2( t)  dt  =  f (i)  +  g (i  +  1)  +  h(i  +  2)  , 
^1 


where  f(i)  is  the  area  under  the  first  i  saw  teeth,  g(i  +  1) 
is  the  shaded  area  shown  in  Figure  4  and  h(i  +  2)  is  the  area 
remaining  under  the  curve. 
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Figure  4.  Finite  Production  Model,  Area  Under  the  Inventory 
Curve  at  Echelon  Two. 


From  Figure  4  it  can  be  shown  that 


Q?  [i-1]  Q?V»  Q?[iM]  [i-1]  ti-2] 
-±  +  - -i—  +  — = - - - 


(3-4) 


Q7  tn-i] [n-i-lj 


(3-5) 


From  Figure  5  the  area  g(i  +  1)  is  found  to  be 

-Q?  ln-1]2  Q?^[i2-1] 

g(i  +  1)  =  —Sip - 


Q?  [i  ln-i]  ]  Q?  {  [i— 1  ]  [iji-X  ] } 

+  ~  x - + - n - • 


r~ 

Qi  Qi 

a  -  [i  jr  -  r  ln'1]] 

b  =  -±  tn-13  -  [i-13 

c  =  0l 

d  —  \J;b 

=  Q1[n-1]  -  Q1  [i—  1} 


Figure  5.  Finite  Production  Model.  Area  Under 
i  +  1st  Saw  Tooth. 


B.  OPTIMAL  POLICY 

The  development  of  an  algorithm  for  determining  the  optimal 
inventory  policy  in  the  finite  production  model  follows  an 
argument  that  is  identical  to  the  one  discussed  in  Section  II. 
Since  (3-9)  is  convex  in  for  fixed  n,  Q£(n)  must  satisfy 


3K(Qirn) 
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or 


Qj(n)  * 


“  ^2 
2tAl  n  ]X 

m 

C[I1+I2[n-l]  - 

I2^tn-21 

1/2 


When  (3-10)  is  substituted  for  in  (3-9)  then 


K(n) 


=  j 

w* 


k: 


KIi  +  1t+k212 


K^2[[n-1]  -  |[n-2]] 


n 


1/2 


where  K?  .  =  2XA.I  .C  . 


(3-10 


(3-1) 


ID  — "i-J 

It  is  apparent  that  (3-11)  is  positive  for  all  values  of  n. 
Further , 


lim 

n+1 

K(n)  =  K^x  +  K21  +  K22X  +  K22^.  . 

and  for  very  large  values  of  nf 

1. 

K(n)^K12n1/'2  >  0  , 

2. 

3K  K12  n 

»»  ZF*  ' 

•5 

32K  „  -K12  <  „ 

<3  • 

3n2  ' 

If  the  assumption  that  n  is  an  integer  is  temporarily 
relaxed  the  slope  of  K(n)  can  be  investigated  using 


+  K?-tl-x] 


If  (3-12)  is  positive  when  n  =  1,  it  is  also  positive  over 
all  n  and  it  follows  that  the  optimal  value  of  n  is  n*  =  1. 
If  (3-12)  is  negative  when  n  *=  1,  then  there  exist  values  of 
n  such  that  K(n)  is  less  than  K(l).  However,  since  it  has 
been  shown  that  (3-12)  is  positive  for  very  large  n,  n*  must 
be  finite  and  satisfy 


be  denoted  by  n,  and  it  follows  that  in  the  range  1  S.  n  £  n 
K(n)  is  convex.  For  values  of  n  >  n  K(n)  is  not  convex,  but 
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(3-12)  is  positive.  Therefore  n*  is  the  optimal  non  integer 
value  of  n  and  (3-11)  is  pseudo-convex. 

Since  (3-11)  is  pseudo-convex  over  all  n,  it  follows  that 
it  is  also  pseudo-convex  on  the  integer  values  of  n.  So  if 
n*  is  not  an  integer,  then  the  optimal  value  of  n  will  cor¬ 
respond  to 

K*  =  min[K(n),  K(n)]  . 

The  optimal  policy  algorithm  can  now  be  written. 

2  *2  A2 

1.  Compute  *22'  an<^  K12  * 

2.  If  Kji  K22  +  then  n*  =  1.  Go  to  step  8. 

2  a2  a  2 

3.  If  K2^  >  K22  +  K12  then  compute  n*  using  (3-13) . 

4.  If  n*  is  an  integer  then  n  optimal  =  n*.  Go  to  step 


5.  If  n*  is  not  an  integer  compute  n  and  n. 

6.  Compute  K(n)  and  K(n)  from  (3-11). 

7.  K*  *  min[K(n),  K(n)]  and  n  optimal  is  the  value  of  n 
corresponding  to  K*. 

8.  Compute  Q*  from  (3-10). 

C.  SPECIAL  CASES  OF  ip 
When  ijj  =  X 

AjX  ^CQ-l  A  2X  I2CQ1 

K<n)  =  of  +  -5 - +  Hq^+  2  ’ 
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lim  K(n) 
n-*» 


AXX  U1+I21CQ1 


(3-14) 


Thus  it  is  seen  that  if  the  production  rate  is  equal  to  the 
demand  rate,  the  optimal  policy  is  to  start  the  production 

line  and  never  let  it  stop. 

As  one  would  expect,  when  ■*  05 

A.X  I^CQ^  A2^ 
lim  K(n)  =  * —  +  — 5 — +  KoT 


I2CQ1ln-l] 


and  it  is  immediately  seen  that  this  model  degenerates  to 
the  EOQ  model  of  Section  XI. 


NUMERICAL  EXAMPLE 
Consider  the  problem  where 

11  =  0.50  ,  C 

12  =  0.20  * 

X  =  25  units/year  ? 

ijj  a*  75  units/year 

2  ~2  A2 

First  compute  K^,  K22  and  K^2 
K21  =  5.0  x  105 

K22  =  2.0  x  105 

Ki2  =  5.0  x  104 


=  $100/unit 


=  $50/unit 


=  $200/unit 


=  1.25  x  10' 


=  6.67  x  10 


=  3.33  x  10 


•  C-vjv;  >if^5 


P?".  WWVJi' 
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u^j  nssy*^ 
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Next  compare  K^  with  K^  +  K^— K^  =  5.0  x  105  >  K^  +  K^2  =  1.0  x  105, 

so  compute  n*  using  (3-13).  The  value  is  3.64  which  is  not  an 
integer.  Therefore  n  =  4  and  n  =  3. 

From  (3-11) 

K (n)  =  $685.00, 

K(n)  =  $694.00, 

and  n  optimal  is  4.  Using  (3-10) 

Q|  =  7.2  units 

and 

Q|  =  28.8  units 

E.  PARAMETER  ANALYSIS 

Figures  11  and  12  illustrate  the  effect  of  production  rate 
on  optimal  n.  If  the  holding  cost  at  echelon  two  is  small 
then  optimal  n  was  relatively  insensitive  to  a  change  in 
production  rate  (Figure  11).  However,  the  results  of  Figure 
12  show  that  this  sensitivity  increases  as  I2  increases. 

The  responses  of  this  model  to  parameter  variations  were 


identical  to  the  responses  observed  in  Section  II. 


IV.  THE  TWO  ECHELON  MODEL  WITH  BACKORDERS  PERMITTED 

A.  MODEL  FORMULATION 

The  EOQ  and  finite  production  models  were  based  on  the 
assumption  that  all  demands  would  be  immediately  satisfied. 
This  assumption  will  now  be  relaxed  in  order  to  investigate 
a  system  in  which  all  demands  are  ultimately  satisfied,  but 
where  it  is  permissible  to  accumulate  backorders  at  the  lower 
echelon.  No  backorders  are  allowed  in  echelon  two. 

In  this  model  demands  which  occur  when  the  lower  echelon 
is  out  of  stock  are  backordered  against  future  procurement. 
When  the  procurement  arrives  these  backordered  demands  are 
met,  and  any  excess  quantity  is  placed  in  stock.  Each  back¬ 
order  at  echelon  one  results  in  a  cost  of  the  form  it  +  it, 
where  u  is  a  fixed  charge  and  it  is  a  time  dependent  cost. 

As  noted,  if  there  were  no  costs  associated  with  incurring 
backorders  [Hadley  and  Whitin] ,  then  it  would  be  optimal  to 
never  have  any  inventory  on  hand.  Conversely,  if  the  back¬ 
orders  are  sufficiently  expensive,  then  the  optimal  policy 
would  be  to  never  incur  any.  However,  for  an  intermediate 
range  of  backorder  costs,  it  will  usually  be  optimal  to  incur 
some  backorders  toward  the  end  of  the  cycle. 

With  the  exceptions  noted  in  the  preceding  paragraphs,  the 
model  is  predicated  on  the  same  assumptions  used  in  Section  II 

Define  s  as  the  number  of  backorders  at  echelon  one  when 

an  order  is  received.  Then  from  Figure  6  it  follows  that  the 
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on  hand  inventory  at  echelon  one  immediately  after  receipt  of 


an  order  is  (Q^  -  s)  units.  It  is  assumed  that  >  s.  Then 
the  echelon  one  on  hand  inventory  varies  from  (Q1  -  s)  to 
zero. 


r~ 


Figure  6.  Two-Echelon  Bc.ckorders  Allowed  Model. 


The  optimal  policy  continues  to  be  that  which  minimizes 
average  annual  variable  system  cost  subject  to  the  model 
constraints.  Since  the  model  does  not  permit  backorders  at 
the  second  echelon,  the  average  annual  variable  cost  at 
echelon  two  will  be  identical  to  equation  (2-4)  ,  which  for 
convenience  is  restated  here. 


(4-1) 
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annual  variable  cost,  K,  as  a  function  of  n  and  the  system 
parameters.  After  regrouping  the  terms  this  expression  can 
be  written 


K(n) 


1/2  1/2 
tlCtljC+TTj-lJc2]  [2AXtI1C+Tr]-Tr2A2]  I-jCttX 


PjC  +  TT] 


I,C+ 


l'“~7T 


(4-10) 


Equations  (4-8) ,  (4-9),  and  (4-10)  are  valid  only  if 
2  ? 

2AXIC  -  it  X  ‘  >.0.  Otherwise  an  analysis  similar  to  that  of 


Hadley  and  Whitin  is  required  for  the  case  when  tt  and  tt  are 
positive.  Such  an  analysis  is  beyond  the  scope  of  this  thesis. 
Because  of  the  complex  form  of  K(n)  given  by  (4-10)  as  a 
function  of  n  no  attempt  was  made  to  evaluate  whether  or  not 


K(n)  was  convex  or  pseudo-convex  in  n  when  tt  and  tt  are  both 
positive.  However,  an  example  is  presented  later  in  which  an 
optimal  value  of  n  >  1  is  obtained  for  tt  and  tt  positive. 

If  tt  =  0  (4-10)  reduces  to 


I1C 


1/2 


K (n)  =  [2AXIC  -  2AXI1C  -—•]  , 


(4-11) 


To  expand  (4-11) ,  let 


Then, 


K  ij  =  2XVjC 


1  2  2  K21X 
K(n)  =  -±j2  Knx  +  ki2x  [n-1]  +  - 


(4-12) 


,  K|2Xtn-13  -2  t  ^  >4  h  c  1/2 

n  11  1  n 
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and  since  x  =  I^C  +  u,  K(n)  reduces  to 


K  (n)  = 


2  2  * 

K11  *  K21  *  2 

-H-r  +  — - -  +  Kl2  [n~l] 


jJE^C+n  nfl^C+i1] 


K^2!n-lJ 


n 


1/2 


(4-13) 


Following  the  analysis  of  Section  II,  it  is  immediately 
seen  that  (4-13)  is  positive  for  all  feasible  values  of  n,  and 
that 


lim  K(n)  = 
n+1 


2  -  2  -  1/2 

KL*  +  K21* 


I,C  +  IT 


Also,  as  n  gets  very  large 

1.  K(n)-*K12  n1/2  >  0  , 


2. 


3K  +  K12  >  n 

^  S?7I  °  ' 


3. 


32K  (n)  ^~K12  ,  n 

— T~  r37?  0  • 


3n 


2n‘ 


Assume,  for  the  sake  of  argument,  that  n  is  continuous 
for  all  n  >,  1.  Then 


32K 


_K2  K2 

*  21  .  „2  ,  *22 

[I^C+tt  ]  n  n' 


3n 


K  (n) 


1/2 


(4-14) 


and  it  is  obvious  that  for  n  =  1  (4-14)  is  non  negative  if  and 
only  if 


36 


ItiitliMiiiiiiiKfc 


i  5 

:  i 


,  * 

>< 


*  'i 
1 


i.  Kiia  WSbJ'A*!  ^St- 


:  vr  '^T.JS.TT 


K21  *  i2 


2  ^1^ 

+  K~o  [1  +  -4-]  . 

TT 


(4-15) 


If  (4-15)  holds,  then  (4-14)  is  positive  for  all  n.  This 
implies  that  over  all  n,  the  slope  of  K(n)  is  positive,  and 
n*  =  1. 

IF  (4-14)  is  negative  at  n  =  1,  then  there  exist  values 
of  n  such  that  K(n)  <  K(l).  However,  it  has  been  shown  that 
for  very  large  n,  (4-14)  is  positive,  and  it  follows  that  n* 
is  finite  and  satisfies 


f  -  ° 


or 


n*  = 


K21:  -  K22  l1!0  + 


IT) 


K12fIlC  + 


TT] 


1/2 


(4-16) 


To  complete  the  argument,  note  that  when  (4-14)  is 
negative  at  n  =  1, 


2 

3  k 


3n 


2 


>  0 


over  the  range  of  n  values  from  1  to  beyond  n*.  Let  the  value 
of  n  where 


2 

3  k 


=  0 


3n 


be  noted  by  n,  and  it  follows  that  K(n)  is  convex  over  the 
range  1  <  n  <.  n.  For  values  of  n  >  n,  (4-13)  is  not  convex, 
but  in  this  range  (4-14)  is  always  positive.  Therefore  n* 
given  by  (4-16)  is  the  optimal  non  integer  value  of  n.  The 
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reader  should  again  note ' that  the  above  arguments  indicate 
that  (4-13)  is  pseudo-convex. 

Since  K(n)  is  pseudo-convex  over  all  n,  it  is  also  pseudo- 
convex  on  the  integer  values  of  n.  Therefore,  if  (4-16)  does 
not  yield  an  integer  value  then  the  optimal  value  of  n  can  be 
found  by  evaluating  K(n)  and  K(n)  as  in  Section  II. 

An  algorithm  for  finding  the  optimal  policy  (in  the  case 
where  tt  =  0)  can  now  be  stated. 

2  2  2 

1.  Compute  K21,  K22'  an<^  K12  * 


I  C 

2.  Test  K21  <  [K^2  +  K22] [1  +  -i~]  . 

IT 

holds  then  n  optimal  is  1.  Go  to  step  8. 


If  the  inequality 


I  C 

3.  If  K21  >  [K?2  +  K22] [1  +  -—-)  then  compute  n*  using 


(4-16) . 


4.  If  n*  is  an  integer  then  n  optimal  =  n*.  Go  to  step  8. 

5.  If  n*  is  not  an  integer  compute  n  and  n. 

6.  Compute  K(n)  and  K(n)  from  (4-13). 

7.  K*  =  min  [K  (in)  and  K(n)]  and  n  optimal  is  the  value  of  " 
n  corresponding  to  K*. 


8.  Compute  Q£  from 


2X  [Al  ♦  5*] 


Q* (n)  *  -- 

C[l1l2C[n-l]  +  Ij*  +  I2tt  [n-l]] 


and  s*  from 


s*(n)  = 


IXC  Q*  "1 1/2 

A 

I.C  +  IT 
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C.  NUMERICAL  EXAMPLE 

Consider  the  problem  where 


*1  - 

0.75 

TT  = 

$5.00 

X2  = 

0.25 

/\ 

IT  = 

$100. 00/year 

A1  " 

$25.00 

X  * 

100  units/year 

a2  =• 

$100.00 

c  = 

$100/unit 

2 

Begin  by  computing  =  2XA^IjC. 


K21  =  1-5  x  10° 

K^2  ~  5.0  x  105 

K^2  =  1.25  x  105 

I  C 

Next  compare  K21  <  [K^  +  K22] [1  +  . 


The  inequality  does 


not  hold  so  compute  n*  using  (4-17).  The  value  is  1.69  which 
is  not  an  integer.  From  (4-13) 


K (n)  «  $735.00 

and 

K(n)  =  $620.00 


Therefore  optimal  n  =  1.  The  values  of  Q£  ,  Q*  ,  and  s* 
are  now  found  to  be 


Q£  =  12.5  units 

Q£  -  12.9  units 

s*  =  8.16  units. 
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D.  PARAMETER  ANALYSIS 


As  noted  previously,  computational  difficulties  preclude 
an  analytical  treatment  of  the  general  form  of  the  backorder 
model.  However,  the  response  of  the  model  to  parameter 
variations  can  be  approximated.  This  can  be  accomplished  by 
considering  separately  the  two  functions  (3-1)  and  (3-2)  and 
ignoring,  for  now,  the  integer  restriction  on  n. 

Equation  (4-1)  can  be  rewritten 

A,X  I-CtCU  -  Q.] 

K„  =  ~~~  +  ^ - —  , 


and  since  the  second  term  is  never  negative,  it  is  obvious 
that  K2  is  convex  in  Q2  .  Using  the  Hessian  it  can  be  shown 
that  (4-2)  is  convex  in  Q.  and  s  if  and  only  if 


2AXX  [1^  +  tt]  >  [irX  ]  2 


(4-18) 


Then 


sr +  - +  kms  +  T"1 


V  X2C[Q2  "  Ql3 


(4-19) 


is  convex  if  and  only  if  (4-18)  holds. 

It  can  be  shown  that  when  (4-19)  is  not  convex  the  optimal 
inventory  policy  will  be  realized  when  backorders  are  not 
allowed  (s*  =  0).  Thus  if  (4-18)  does  not  hold,  the  problem 
can  be  solved  by  setting  s*  =  0  and  using  the  algorithm 
described  in  Section  II  to  get  Q£  and  n*. 


For  the  case  where  (4-19)  is  convex,  optimal  Q^,  Q2  and 
s  can  be  found  from 


3K 


3K 


3K 

Ts 


=  0  . 


(4-20) 


Equation  (4-20)  will  be  satisfied  when 


Q1  " 


2\  [nA^  +  A2)  [I^C  +  tt]  —  nir 


2X  2 


^n[I.jC  +  it  +  I^C  [n-1]  +  I2CiT[n-l] 


1/2 


(4-21) 


s* (n)  * 


and 


I jCQ*(n)  -  ttX 


(4-22) 


IjC  +  tt 


Q^(n)  =  nQj(n)  , 


(4-23) 


and  the  optimal  policy  can  be  found  from  the  following 
algorithm. 

1.  Set  n  =  1  . 

2.  Compute  Qj(n),  s*(n),  Q^(n)  and  K(n)  using  equations 
(4-21)  through  (4-23)  and  (4-19). 

3.  Set  n  =  2  and  repeat  step  2. 

4.  If  K(l)  <  K(2),  stop;  n*  =  1  and  the  optimal  policy 
is  known. 


5.  If  K(2)  <  K (1) ,  set  n  =  3  and  repeat  step  2. 

6.  Continue  solving  K(n)  by  increasing  n  in  steps  of  1 
until  K(n+1)  >,  K(n).  Stop  as  soon  as  K(n)  <  K(n+1);  n*  = 
current  val  le  of  n  and  the  optimal  solution  has  been  found. 
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',j^'’rr5vjf=tyv^w’5’jyv^  . 


'■?,  •  '  "*v**?.-.  ,-v^-j? 


''•  v*-  *>.sf,STI 


The  algorithm  described  in  the  preceding  paragraph  was 
used  to  observe  the  effect  of  parameter  variations  on  optimal 
n.  In  general,  the  results  were  similar  to  those  observed  in 
the  EOQ  model.  For  example,  n  optimal  varied  inversely  with 
A1  and  I2,  but  directly  with  A2  and  1^. 

The  sensitivity  of  n  to  the  various  parameters  was  highly 
dependent  on  the  selected  values  of  the  backorder  costs.  For 
example,  Figure  13  illustrates  the  relationship  between  optimal 
n  and  I2  in  a  situation  where  the  backorder  cost  is  small. 

Figure  14  illustrates  this  same  example  with  a  high  time  de¬ 
pendent  backorder  cost.  The  difference  in  the  sensitivity  is 
obvious . 

Figures  15  and  16  illustrate  the  relationship  of  optimal  n 

to  a  change  in  A2>  In  Figure  15  the  backorder  costs  were  low, 

and  optimal  n  was  insensitive  to  a  change  in  A2>  In  this 

example  optimal  n  remained  at  two,  and  s*  increased  as  A2 

increased.  When  the  time  dependent  backorder  cost  was  high 

optimal  n  varied  directly  with  A2  (Figure  16) .  However,  the 
*  * 

change  in  optimal  n  was  not  nearly  as  large  as  in  the  EOQ  model. 

Figures  17  through  20  illustrate  the  relative  magnitude  of 
QJ  ,  Q|  ,  and  s*  as  A2  and  I2  are  allowed  to  vary.  The  breaks 
in  these  curves  occur  at  points  where  optimal  n  changes. 


ft 
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V.  RECOMMENDATIONS  FOR  Fl|T*  IER  STUDY 

The  analysis  of  the  backorders  permitted  model  should  be 
completed  for  all  values  of  the  backorder  costs  it  and  tt  in¬ 
cluding  the  case  where  ir  >  0  and  n  =  0.  Hadley  and  Whitin 

discuss  this  case  in  their  single  echelon  development  and 
conclude  that  if  tt  =  0  then  s*  is  either  0  or  infinite. 

Initial  investigations  indicated  that  in  the  two  echelon  model 

A 

there  could  exist  a  finite  value  of  s*  >  0  when  it  =  0.  How¬ 

ever,  this  investigation  was  not  completed  and  no  conclusions 
were  reached. 

It  is  felt  that  the  analytical  argument  used  throughout 
the  paper  could  be  applied  to  the  general  backorder  case  where 
tt  and  tt  are  ooth  greater  than  zero.  However,  the  complexity  of 
the  equations  would  greatly  complicate  this  development. 

The  models  discussed  in  this  thesis  should  be  extended  to 
more  than  one  activity  at  each  echelon.  If  it  can  be  assumed 
that  all  activities  within  a  given  echelon  order  at  the  same 
time,  then  it  is  particularly  easy  to  include  multiple  activ¬ 
ities  within  an  echelon.  For  example,  assume  that  there  are 
K  activities  in  echelon  one.  Then  in  the  EOQ  model  would 
take  the  form 

TliCQU 
2 

If  the  values  of  the  parameters  A^,  X^,  1^,  and  C  are  constant 
for  all  i,  i  =  l,...,k,  then  the  expression  for  reduces  to 
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the  trivial  case  of 


K«  =  k  Kj.  , 

where  is  found  in  Section  II.  The  expression  for  K2 
remains  unchanged. 

The  problem  of  multiple  activities  at  the  top  echelon  is 

i 

usually  not  a  very  interesting  one,  since  each  top  echelon 
activity  is  usually  responsible  for  supplying  specified  activ¬ 
ities  at  the  next  lower  echelon.  Therefore,  unless  the  model 
provides  for  lateral  resupply  actions,  the  problem  of  multiple 
activities  at  the  top  echelon  can  be  reduced  to  the  sum  of  a 
series  of  independent  problems,  each  containing  one  top 
echelon  activity. 

A  related  problem  which  should  be  investigated  is  the  ex¬ 
tension  of  the  models  to  more  than  2  echelons.  In  the  case 
where  each  echelon  is  limited  to  one  activity,  this  extension 
would  not  be  difficult.  Indeed,  given  the  assumption  that  there 
is  an  integer  ratio  between  the  quantity  ordered  by  successive- 
echelons,  the  model  would  probably  take  the  form 

K  =  J  K.  , 

j=l  3 


where  the  relationship  between  the  Qj's  is 


Q2  =  ni  *  Qi 


Q3  =  n2  •  Q2  =  n2  •  n1  •  Qx 


Q .  —  n .  ,  •  Q.  .  =  II  n. Q.  , 

3  3i  3  l  i=1  i  1 

and  nK  is  an  integer  for  all  k,  k  =  l,...,j-l.  Because  of  the 
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integer  property  of  the  n^'s,  the  solution  technique  outlined 
in  Chapter  II  could  be  used  to  solve  recursively  for  K. 

If  the  m-echelon  problem  allows  k^  activities  at  the  ic 
echelon,  i  =  l,...,m,  then  for  large  values  of  k^  and  m,  the 
sheer  magnitude  of  the  problem  would  make  its  solution  ex¬ 
tremely  difficult.  In  fact,  it  is  doubtful  that  the  proce¬ 
dures  suggested  in  this  paper  could  be  utilized  for  a  problem 
of  this  type.  However,  if  the  model  is  small  enough,  one 
could  at  least  get  a  "feel*1  of  the  model's  behavior. 


VI.  SUMMARY 

Three  of  the  deterministic  models  of  Hadley  and  Whitin 
have  been  restructured  as  two  echelon  models.  Equations  for 
the  average  annual  variable  cost  are  derived,  and  standard 
mathematical  programming  techniques  are  utilized  to  find  the 
optimal  inventory  policy.  The  optimal  policy  is  defined  as 
that  which  minimizes  average  annual  variable  cost  subject  to 
the  constraints  of  the  model. 

General  solutions  are  obtained  for  the  EOQ  and  finite 
production  models.  Because  of  the  complexity  of  the  cost  equa^ 
tion  the  analytical  solution  of  the  backorders  permitted  model 
is  limited  to  the  case  where  tt  =  0.  However,  a  technique  is 
developed  which  can  be  used  to  find  the  optimal  policy  in  the 
general  case. 

A  parametric  analysis  is  conducted  for  each  model  in  order 
to  study  the  behavior  of  the  models  under  parameter  variations 
The  behavior  of  the  curves  of  n  optimal  for  the  EOQ  and  finite 
production  models  were  identical  except  that  the  finite  produc 
tion  model  enjoyed  a  significant  cost  advantage.  The  behavior 
of  the  backorder  model  was  very  dependent  on  the  magnitude  of 
the  backorder  cost.  For  example,  if  the  model  had  a  low  time 
dependent  backorder  cost  then  parameter  variations  had  little 
effect  on  the  optimal  behavior  of  the  model.  However,  if  the 
backorder  cost  was  set  sufficiently  high  then  the  model's 
behavior  was  similar  to  that  of  the  EOQ  model.  Figures  7 
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ijW/t. 
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Figure  8 .  Optimal 


Figure  10.  Optimal  n  :  I2  •  EOQ  Model 


Optimal  n  :  I roduction  Rate.  Finite  Production  Model 


Figure  12.  Optimal  n  :  Production  Rate.  Finite  Production  Model 


Figure  13.  Optimal  n  :  I2  with  Low  Backorder  Cost.  Backorders  Allowed  Model 


Optimal  n  :  I-  with  High  Backorder  Cost.  Backorders  Allowed  Model 


Figure  15.  Optimal  n  :  A,  with  Low  Backorder  Cost.  Backorder  Allowed  Model 


Figure  16.  Optimal  n  :  h~  with  High  Backorder  Cost.  Backorder  Allowed  Model 


Backorder  Costs. 
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Figure  20.  Q£  ,  q*  and  s*  .  j  Backorder  Model  with  High  Backorder  Costs 
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